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ABSTRACT 


The  present  paper  extends  the  method  of  Cheng  and  Trumpler  (s)  to  study 
stability  of  plain  cylindrical  gas  Journal  bearings  of  finite  length.  Both 
equilibrium  and  stability  results  have  been  obtained. 


INTRODUCTION 
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It  has  been  shovm  recently  by  many  authors,  [l  -  loj  that  one  of  the 
most  Important  considerations  In  designing  a  high  speed  gas-lubricated  self¬ 
acting  Journal  bearing  Is  the  Instability  of  the  Journal  under  a  given  operat¬ 
ing  condition.  Intensive  research  In  this  direction  has  led  to  a  number  of 
significant  contributions  In  the  past  five  years. 

The  first  attempt  on  this  stability  problem  was  accomplished  by  Sternllcht 
and  Rentzepls[2  ].  Their  analysis  Includes  only  partial  effect  of  the  time 
variation  of  pressure  in  the  Reynolds  equation.  The  results  predicted, 
for  the  first  time,  the  existence  of  the  worst  clearance  ratio  at  which  the 
threshold  speed  of  Instability  Is  minimum.  The  same  results  also  Indicate 
the  Stability  Is  enhanced  by  a  large  compressibility  number. 

f6}  (8]  .  . 

Recently  Castelll  and  Elrod,  Ausman,  and  Cheng  and  Trumpler|^3  Jhave 

A  A 

provided  stability  solutions  for  the  Infinitely  long  bearing.  The  results 
of  these  Investigations  show  reasonably  close  agreement. 

A  detailed  summary  and  comparison  of  experimental  data  and  various 
solutions  to  this  problem  are  Included  In  a  recent  work  by  Pan  and  Sternllcht[l^. 
They  also  give  the  results  of  a  stability  analysis  using  the  quasl-statlc 
linearized  ph  stability  solution  of  the  Reynolds  equation  for  bearings  of 
finite  length. 

The  major  obstacle  In  this  problem  lies  In  the  difficulty  to  obtain  an 
accurate  solution  of  the  non-llnear  Reynolds  equation  with  the  time  dependent 
term.  The  method  of  B.  G.  Galerkln  applied  to  the  function  ph  ( 5  }  has  proven 
to  be  a  very  effective  way  to  handle  equations  of  this  kind.  It  has  the 
main  advantage  of  reducing  the  Rejmolds  equation  directly  from  the  partial 
differential  equation  to  a  system  of  first  order  ordinary  differential  equations 
which  together  with  the  equations  of  motion  of  the  Journal  yield  quite  readily 
to  stability  analysis  of  the  system.  It  also  can  be  applied  to  bearings  of  more 
complicated  configurations  such  as  partial  arc  bearings,  so  long  as  a  dependable 
quantitative  description  of  the  ph  function  can  be  guessed.  . 

The  present  report  Is  essentially  the  extension  of  the  work[s]to  the 
stability  study  of  finite  Journal  bearings  using  Galerkln's  method. 
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I,  GOVERNING  EQUATIONS 

Considering  a  finite  length  journal  and  bearing  system  as  shown  in 
Fig,  1,  the  equation  governing  the  pressure  in  the  gas  film  is  the  iso¬ 
thermal  Reynolds  equation. 

5{  l|>  *  Is  H)  ■  A  <1  -  ^  ^  .  <1) 
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Assuming  the  journal  is  perfectly  balanced  and  that  it  rotates  at  a  con¬ 
stant  rotational  speed,  the  equations  governing  the  translational  motion 
(as  opposed  to  the  conical  motion)  of  the  journal  become 
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Let  \|r  -  PH,  equation  (1)  becomes 


-  3 


(4) 


In  order  to  apply  the  method  of  Galerkln  to  equation  (4),  we  assume 
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where  f  >  JJ  t 

Using  Galerkln'a  method  for  m  -  2  and  n  ■  1,  equations  (2),  (3)  and  (4) 
are  reduced  to 
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The  detail  derivations  of  equations  (6)  to  (12)  are  Included  in 
Appendix  I. 
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EQUILIBRIUM  SOLUTION 
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By  setting  all  derivatives  equal  to  zero  in  equations  (6)  through 
(12) ,  we  obtain  readily  the  equations  governing  the  static  equilibrium 
position  of  the  Journal.  Numerical  solutions  of  equation  (8)  through 
(12)  for  obtained,  using  Newton-Raphson' s 

method  [itJ  .  Once  Cj^^.A^^^^etc,  are  determined,  the  expression  for  the 
load  and  attitude  angle  can  be  derived  from  eqs,  (6)  and  (7).  They  are 


I  ks*  ^  Kb»  5 tea 

cko  -  -  .  .  (1^) 

(k'Ba-  Kso) 

The  results  of  the  equilibrium  solution  are  tabulated  in  Tables  1  to  4 
and  plotted  as  the  load  and  attitude  angle  charts  in  Figs.  2  to  9. 

Comparison  between  the  present  solution  with  two  finite  difference 
numerical  solutions  [.1^  [l2]  Is  also  shown  in  Figs.  lO  to  13. 

In  general,  the  present  results  compare  quite  favorably  with  finite 
difference  computer  results.  The  largest  error  Introduced  by  using 
Galerkln's  approximation  appears  to  fall  in  two  regions: 

1.  High  €  and  low  A  region:-  In  this  region,  the  present  theory 

predicts  a  lower  load  and  a  higher  attitude  angle  than  finite 

difference  solutions.  This  type  of  error  Is  due  to  the  trunca¬ 
tion  error  of  the  assumed  function  in  the  0  direction.  The 

same  characteristic  error  was  also  observed  In  the  analysis  of 

infinitely  long  bearings. 

2.  Large  ^  and  high  A  region;-  In  this  region,  the  Galerkin  method 
gives  a  slightly  lower  load  and  attitude  angle  comparing  to 
finite  difference  solutions.  The  reason  for  this  error  is  that 
only  one  cosine  term  Is  assumed  In  the  axial  direction.  It  is 


conceivable  that  the  results  will  be  much  Improved  l£  more 
terms  are  considered  In  the  axial  direction.  However,  the 
algebraic  complexity  Introduced  by  assuming  more  terms  In  the 
axial  direction  would  also  become  much  more  formidable. 


IV. 


STABILITY  OF  THE  EQUILIBRIUM  SOLUTION 


The  stability  of  the  equilibrium  solution  can  be  Investigated  by 
considering  a  small  perturbed  displacement  of  the  journal  center  and 
small  variations  of  the  pressure  coefficients,  C,  A,  etc.,  which  are 
defined  by 
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Substituting  (15)  Into  equations  (6)  to  (12),  and  neglecting  all  non¬ 
linear  terms,  we  obtain  the  linearized  equations  of  the  dynamical  system 
which  can  be  written  In  the  matrix  form  as 


AX  +  I 


dr 


-  0 


(16) 


where  the  column  matrix  X  represents 


'll. 

^12. 

*12 


(17) 
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I  is  the  unit  matrix  and  A  Is  a  9  x  9  coefficient  matrix  whose  elements 
are  expressed  In  terms  of  the  parameters  P^,  A,  and  equilibrium 
solutions  e^,  Cio*'^llo  ^  cnniplete  list  of  these  coefficients 

Is  given  In  Appendix  II. 

The  characteristic  equation  for  the  system  represented  by  equation 
(16)  Is  a  nlneth  degree  polynomial  whose  coefficients  can  be  determined 
numerically  by  using  Danielewsky ' s  method  ^is]  .  After  the  characteris¬ 
tic  equation  Is  determined,  the  stability  of  the  equilibrium  solution 
can  be  Investigated  by  applying  Routh's  Criterion  ^14^  to  the  coeffic¬ 
ients  of  the  polynomial.  For  a  given  set  of  parameters  and  A, 

It  Is  found  that  there  exists  a  threshold  speed  below  which  the  Journal 
position  Is  stable  and  above  which, unstable.  The  numerical  calculation 
of  the  polynomial  coefficients  as  well  as  the  threshold  speeds  are  per¬ 
formed  on  an  IBM  7090  digital  computer.  The  results  of  the  threshold 
speed  for  different  values  of  and  A  are  tabulated  In  tables  1  to  4 

and  also  are  plotted  as  the  stability  charts  In  Figs.  14  to  17. 

The  whirl  frequency  of  the  system  at  the  threshold  of  stability  can 
be  determined  by  extracting  the  roots  of  the  characteristic  polynomial 
at  the  threshold  speed.  The  complex  conjugate  roots  which  have  a  nearly 
zero  real  part  are  Identified  as  the  whirl  frequency,  since  its  real  part 
Is  at  the  threshold  of  becoming  positive.  Complete  whirl  frequency  maps 
are  calculated  for  ^  equal  to  1  and  1/4  and  are  shown  in  Figs.  18  and  19. 

The  parameters  used  In  the  analysis  are  convenient  for  computation. 

They  cannot  be  directly  used  for  design  purpose,  since  both  the  parameters 
/  MCN*^ 

/and  A  Involve  the  rotational  speed.  Another  set  of  stability  para¬ 
meters  more  suitable  for  designing  purpose  was  first  given  by  Rentzepls  and 
Sternllcht.  The  advantage  of  that  plot  Is  that  the  speed  can  be  varied 
Independently  without  changing  other  parameters.  The  conversion  of  the 
present  stability  data  to  C  and  cu  plot  are  shown  In  Fig  20  and  21  for 

^  «  2  and  1/2. 
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V.  DISCITSSICTI  OF  RESULTS 

1.  Equilibrium  Solutions 

The  overall  agreement  between  the  equilibrium  solutions  by  the 
Galerkin  method  and  the  finite  difference  method  is  quite  encouraging, 
considering  the  assumed  functions  for  ph  only  Include  a  few  leading 
terms  of  the  Fourier  series.  The  error  contributed  by  the  Insufficient 
terms  in  the  circumferential  direction  becomes  more  pronounced  for  small 
large  e  and  small  A  The  error  of  load  factor  for  ^  =  1/2,  e  •  0.  8 
varies  from  17.  to  10%  in  the  high  A  region  and  rises  to  approximately 
35%  in  the  low  A  region.  The  error  disappears  rapidly  with  decreasing 
e.  The  second  type  of  error  is  attributed  to  the  insufficient  terms 
in  the  axial  direction.  Since  only  a  single  cosine  term  is  chosen  in  this 
axial  direction,  it  is  expected  that  the  error  will  be  amplified  when 
both  ^  and  A  become  large,  regardless  of  the  value  of  e.  For  ^  »  2, 

€  ■  0. 1,  A  “  100,  this  type  of  error  increases  to  18%  for  the  load 
factor. 

2.  Stability  Curves 

In  general,  the  stability  curves  show  the  same  trend  as  those 
obtained  earlier  for  the  infinitely  long  bearing.  Plotting  the  threshold 
speeds  for  small  values  of  A  (on  regular-scale  graph  paper),  one  can  show 
that  the  threshold  speed  tends  to  approach  zero  as  A  goes  to  zero,  regard¬ 
less  of  the  value  of  e  and  This  fact  verifies  the  theory  that  a  non- 
cavltated  full  film  Incompressible  journal  bearing  is  unstable. 

For  constant  values  of  ^  and  c,  the  stability  curve  rises  rapidly 
with  Increasing  A,  forms  a  hump  in  the  middle  range  of  A  and  finally 
approaches  an  asymptote  for  extremely  high  A  The  location  of  the  hump 
and  the  maximum  threshold  speed  depend  greatly  on  ^  and  e.  However,  the 
asymptotic  speed  at  high  A  for  each  e  seems  to  change  little  with  diffe¬ 
rent  values  of  — . 

D 

For  ^  equal  to  1/2  and  1/4,  s  equal  to  0.8  and  0.9,  and  high  values 
of  A,  the  present  theory  falls  to  predict  a  meaningful  threshold  speed. 

In  this  region,  Routh's  Criterion  gives  a  very  low  threshold  speed  with 


the  corresponding  whirl  frequency  many  times  higher  than  1/2.  It 
is  extremely  difficult  to  determine  whether  this  low  threshold  speed 
actually  exists  in  reality  or  is  caused  by  an  extraneous  root  due  to 
the  truncation  error  in  the  assumed  function.  This  region  represents 
either  extremely  heavily  loaded  or  low  ambient  cases  and  thus  is  of 
little  practical  interest.  Such  uncertain  results  in  this  region  are 
discarded, 

3,  Whirl  Frequency  Curves 

The  whirl  frequency  curves  show  that  for  a  fixed  ^  and  e,  the  fre¬ 
quency  approaches  an  asymptote  at  small  A,  reaches  a  minimum  value  In 
the  middle  range  of  A  and  finally  approaches  half-frequency  at  extremely 
high  A.  It  should  be  of  interest  to  compare  the  asymptotic  frequency 
at  small  A  with  that  of  an  incompressible  experiment. 

4,  Comparison  with  an  Earlier  Theory 

The  comparison  of  the  present  stability  curves  and  Pan  and  Stern- 
licht's  Quasi-Static  theory  {16}  shows  general  agreement  in  trend  for 
large  values  of  A  (see  Figs.  22,  23).  On  the  C  and  cu  plot,  they  seem 
to  give  similar  slopes  for  a  constant  P^,  but  the  present  theory  predicts 
a  higher  threshold  speed  in  this  region.  At  small  A  region,  the  two 
theories  give  quite  different  slopes  with  the  present  theory  predicting 
a  more  conservative  threshold  speed. 

5,  Comparison  with  Experiment 

Unfortunately,  there  exist  only  a  few  experimental  data  for 
comparison.  Fig.  22  shows  some  comparison  with  results  made  by  Stern- 
llcht  and  Winn  r9j,  and  Whit ley -Bowhill-McEwan  f?}.  For  ^  *  1  and  P  = 
0.2,  agreement  is  excellent  in  small  A  region,  but  less  satisfactory 
in  high  A  region.  For  ^  “  2  and  P^  »  0.1,  agreement  is  very  good  with 
Sternllcht-Wlnn' s  upper  curve  (the  upper  curve  is  the  spontaneous  in¬ 
stability,  while  the  lower  data  is  obtained  by  giving  an  impulsive  force 
at  successive  speeds,  and  the  speed  at  \:’.^ch  the  stability  cannot  be  re¬ 
stored  is  considered  as  the  thresh  Id  speed).  It  is  surprising  that  the 
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few  date  obtained  by  Whltley-Bowhill  and  McEwan  fall  extremely  close  to 
the  present  stability  curves.  Such  close  agreement  Is  probably  acclden* 
tal. 

Cio) 

The  comparison  with  Reynolds  and  £ross  data  Is  shown  In  Fig.  23. 
These  curves  seem  to  show  a  consistent  trend,  that  the  agreement  Is  very 
close  at  A  values  smaller  than  1,  but  less  satisfactory  In  higher  A 
region. 
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VI.  CCTtCLUSIONS 

1.  The  method  of  Galerkln  provides  an  effective  technique  to  solve  the 
time-dependent  non-linear  Reynolds  equation,  and  It  lends  Itself  very 
conveniently  In  determining  the  stability  curves  for  various  operating 
parameters. 

2.  The  stability  results  using  Galerkln' s  method  show  the  same  trend 
Indicated  by  earlier  work  (see  Fig.  22).  For  example,  a  minimum 
threshold  speed  Is  shown  to  exist  In  all  theories  when  varying  the 
clearance  under  constant  load.  Also,  the  effect  of  ^  Is  similar 
as  given  by  Fan  and  Stemlicht. 

3.  Comparison  of  the  present  stability  results  with  available  experi¬ 
mental  data  Is  encouraging  (see  Figs.  23,  24). 
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VII.  RECroiMENDATIONS 

1.  More  experimental  data  Is  needed  to  establish  a  consistent  cor* 
relation  with  theoretical  results  such  as  given  In  this  paper. 

2.  The  non-llnear  Galerkln  method  should  be  applied  to  analyze  the 
partial  arc  bearings  to  obtain  both  equilibrium  and  stability 
results.  (This  analysis  Is  presently  under  way  at  Mil.) 

3.  The  response  of  the  cylindrical  gas  Journal  under  finite  dynamic 
load  (  particularly  Impulsive  load)  can  be  studied  by  performing 
tlme*wlse  Integration  of  the  simultaneous,  ordinary  differential 
equations  on  either  an  analog  or  digital  computer. 

4.  The  present  analysis  can  be  readily  extended  to  consider  the  stabi¬ 
lity  problem  for  the  conical  mode  of  a  symmetrical  rotor  bearing 
system.  It  also  can  be  extended  to  study  the  dynamics  of  a  non- 
symmetrlcal  rotor.  (This  analysis  Is  presently  under  way  at  MTI. ) 
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APPENDIX  I 


Derivation  of  eg.  (8)  to  eg.  (12) 


•4.  a 


Substituting  eq.  (15)  Into  the  left  side  of  eq.  (4),  we  have 
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[(Jl*>^-Oo>rs  z(£4n-0l  ]  J 

V  M  ^  r 

~~  ^  ^  ^  4  ^  fm40(m-0  Six  (m-hOQ  4  -  (m*0(r>'i-OSii^(m-i)B 

jt-l  (Hi,  ft-.,  I 

[coSz(^4n.i)S  4  J  t  (Zft-i)S  ( >  +  <=  S.  U  hn  $ 

[(jt4rt-i)Cffiz(^4n-i)f  i-  (^-n)c^2(i-4>)i  ]j 
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H  N 


■4'!  '»*'  7  ^  ■) 

-  (2n-i)S  (/-f^cuia)  c^AefC^*  i^-i)c>nzfJ*»-i)S  •^U-*i)Ci/:2f£->il^r 

^■-,  Jlz,  n-.,  (  ^ 

-(zn-i)^0*€ccia)  c^i 2.(£-^ r)-i)i4.(J£-'r> )on. 

V  M  A/  J- 

hm\  ^ 

lf»;  /M.-/  «»/ 

■*'2(^-0(»''*-^*l)Cfi(nit4i-i)d  j(rr>-i)(fn-J^*l')C(C(  >)e 

■i-  ^  {f>^*A-  OoJi  (n>*A—0^  ^  (/n* I )( i^-A-  0  ^  )a] 

f^C*Szf£i-n-i)f  -h  £nZI'/-n)iJ-t-  (2n-oS  O 6  ana)  c/rS^S  C^md 
[(i*n-i)c4i  z(i*fi-i)  %  4-  (j£-n  )C<jS  Z(J£-«>)  f  7  j* 

M  H  M  H 

-1222  6«A.  ^  )5—  M(or)--^)  S>r>  ('m-j4.')d 

-i.-l  i-i  m-i  "*^1  -  ^ 

*  2  C'^-i')(>r>-tA*  l)  S-'n  (n)44.+<)$  -  ^  ,)si'r('>^-A*>)  ^ 

■^%(>^+l)(nr)i'^-i')s',n(m-^A—  l)a  -^(rn-*-d(>^--^-i)  Z  (r^-A  ~  0^^ 

[c4nz.(£-*rf-i)f  ^  J  ^  (2>^‘i)S'  (i-i-  oaikr>9 

ld*n->)  CinifJl-^i^-Or  ■>■(£•  1^)0^  2(£-n)f  J  j 

M  M  M  rj  r 

"XZZ^  S'n/'rr>A)a  -^^^A)rr, 

■^>.  I  ni»i  n»  I  ( 

+  2  (iv-iym-yi-t-  /')  J-'rt  )$  *  ^  ( m-i)('VlfAi‘')  s<n  (<>^■*>4*/)  9 

*  ^  (m*'  )('^'^'4-l)  S'^  ("^*Ai-i)  &  -t  2 

lc4S2CJ*n-i)f  z  Cjt\  Z(l-r>  )f  J  4  (Z'^-l)  S\/i-&c^e)  onA^  5.'< 

f  (i-n)  (VI  ]  j" 


/fl  )n0 
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ti  ti 

WJ»/  /‘/  nti*  •'‘I  ^  . 

+  I  f /»-')(«>  ■'■■4  ♦()<>s/''«-Ni+/)d  -*  l('^~0(n~’y)-l)  ^ 

■^^(^rfi4i')(^*m-  t)  Cin(4im-i)Q  -t-  >^+')  f4->^*>)fij 

[  osM^{jt-^n~i')  f  f  c^  z(^->^)i  ]  *  (Zi^-i)  S  ( i-t-e  ccne)S"r,^d s>nmfi 
2rAK7-/)f  ^  ]  j  ^  A1  ) 


Letting 


then  the  Galerkln  procedure  gives  the  following  set  of  equations: 


•*  J  ?->r 


J  I  (P(f^  <UrS(z/l-l)i  Jl  de  =  O 


/I  ..  /^  ■  -  ■  A/  (  A3  ) 


2-n- 


y  y  e^l>B  di <^9  *  0 


A.~  I,"  ■■  d 
f>=  h  -  -■ 


i'll  t-tr 


j  J  Of  If) 


CfsCz/i-i)^  Si'n^S  df  d9  ~  o 


A.-^  I, - rJ 


Substituting  (Al)  Into  (A3),  (A4) ,  and  (A5),  and  letting  M  »  2  and  N  =  1, 
we  have  eq.  (8)  to  eq.  (12). 


i 

4 
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Derivation  of  eg.  (6)  and  eg.  (7) 

After  substituting  \|r  Into  the  fluid  force  terms 
In  the  equations  of  motion,  the  radial  force  becomes 


i- 5, t  S‘»Z6)  J6 

s.  I^^C,  *  Afj  All*  ( Ki-  )  /4,j,  J 


(  A6  ) 


Similarly,  the  tangential  force  becomes 


Fs 


+  7.d^  dd 


KTs  Bn  ) 


(  A7  ) 


where 


Z'TT 


inr_ 

S*  ! 


K6  =- 


t'irt 

S(S*i) 


^6 

s 


znr6 

(S*  O'- 


-  (/-  €'^') 


(  A8  ) 
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APPENDIX  II 


LINEARIZATION 


Equations  (6)  to  (12)  are  linearized  by  expanding  each  term  In  Taylor's 
series  with  respect  to  the  equilibrium  solution  and  only  keeping  the  linear 
term.  The  procedure  Is  straightforward;  therefore,  only  two  examples  will 
be  given  here  and  the  rest  results  are  listed  as  the  coefficients  in  table  5 

1.  Linearization  of  F 

r 


S'  C,,  V-  V.  C <(.,  -/Ta.  )  /4/z. 

*  { H'  ■="  * 

■+  C,  f  Afi  +  (  *  KBo  ) 


(  A9  ) 


2.  Linearization  of 


An  A/t* 

€  (^AnAn*  Bn  )  ^  ^  ^  C  ^Ho  A?.  ’*■  ^Ho  ^/lo  ) 

^(A//oA,2„*  S/nB/z,}  ^BnA>  ) 


(  Aio  ) 


where  the  subscript  o  represents  the  equilibrium  solution. 
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NCMENCLATURE 


A 

C 

C 

D 

e 

F 

F 

r 

F 

h 

H 

I 

m,  n,  k,  1 
M,  N 
P.  q.  r 
L 

•^3 

K, 


Coefficient  matrix  of  the  stability  equation. 

Pressure  coefficients. 

Radial  clearance,  In.  (also  used  as  a  pressure  coeff. }. 
1/5  „ 


P  2 

/I®)  /_*) 
If  ■'  ^6u'' 


—  dimensionless  clearance. 

K 


Diameter  of  the  Journal,  In. 

Eccentricity,  In. 

Load,  lb. 

Radial  load,  lb. 

Tangential  load,  lb. 

Film  thickness.  In. 

Dimensionless  film  thickness. 

Unit  matrix. 

Index  for  the  pressure  coefficients. 

Index  for  the  order  of  approximation  of  ph  function. 
See  equations  (A3),  (A4)  and  (A5). 

Length  of  the  journal.  In. 

2n 

S(S+1) 

2ff 

S  +  1 

.  ■ 

S(S+1) 


-  21  - 

1 

I 

•^8 

-  -22S-* 

Ts+ip 

1 

1 

1 

5 

n  D 

2  L 

1  • 

M 

Hass  of  the  Journal,  slug. 

P 

Pressure,  psla 

P« 

Ambient  pressure,  psla 

P 

P/Pfl 

F 

P 

m 

p  LD 

B 

R 

Radius  of  the  journal.  In. 

S 

(1  - 

t 

Time,  sec. 

X 

See  definition  (17) 

Z 

Coordinate  In  the  axial  direction.  In. 

e 

Eccentricity  ratio. 

2 

Viscosity,  lb. sec, /in. 

T 

cut 

i  • 

a 

Attitude  angle. 

9 

Coordinates  In  the  circumferential  direction. 

1 

z/R 

8 1  6 

Rotational  speed  of  the  Journal. 

[<a,^  <^1 

^  a  ^  0) 

♦ 

PH 

A 

,R»2 

Pa  ^C^ 

\ 

whirl  frequency 

CO 

• 

n 

P„LD 

a 

i 
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Load,  Attitude  Angle  end  Threshold  Speed  for 


2 


A 

P 

m 

a 

0 

MC  ^ 
oiCy) 

0.1 

0.00789 

86.38 

0.3203 

0.5 

0.03774 

72.50 

0.6356 

1.0 

0.0672 

57.87 

0.7846 

2.0 

0.0997 

38.70 

0.8765 

5.0 

0.1227 

17.83 

0.8697 

10.0 

0.1275 

9.141 

0.8221 

20.0 

0.1288 

4.600 

0. 7900 

100.0 

0.1292 

0.922 

0.7838 

0.1 

0.01596 

86.16 

0.4453 

0.5 

0.07622 

71.55 

0.9115 

1.0 

0.1357 

56.69 

1.1038 

2.0 

0.2033 

37.92 

1.2331 

5.0 

0.2548 

17.56 

1.2235 

10.0 

0.2660 

9.011 

1.1566 

20.0 

0.2691 

4.535 

1.1114 

100.0 

0.2700 

0.909 

1.1027 

0.1 

0.0332 

85.23 

0.6797 

0.5 

0.1593 

67.79 

1.3647 

1.0 

0.2863 

52.03 

1.6046 

2.0 

0.4431 

34.72 

1.7487 

5.0 

0.5916 

16.48 

1.7419 

10.0 

0.6298 

8.489 

1.6262 

20.0 

0.6405 

4.275 

1.5563 

100.0 

0.6440 

0.8568 

1.5381 

0.1 

0.0529 

83.61 

0.9609 

0.5 

0.2620 

61.44 

1.8993 

1.0 

0.4916 

44.41 

2.0997 

2.0 

0.7953 

29.10 

2.1407 

5.0 

1.1491 

14.46 

2.0988 

10.0 

1.2655 

7.5402 

1.9431 

20.0 

1.2997 

3.8035 

1.8444 

100.0 

1.3110 

0.7626 

1.8084 

0.1 

0.0777 

80.56 

1.4765 

0.5 

0.4319 

50.84 

2.901 

1.0 

0.9081 

33.01 

3.0202 

2.0 

1.5684 

20.54 

2.5896 

5.0 

2.4178 

10.80 

2.3519 

10.0 

2.7931 

5.805 

2.1038 

20.0 

2.9146 

2-943 

1.9641 

100.0 

2.955 

0.591 

1.8950 

0.1 

0.0967 

76.81 

2.0391 

0.5 

0.6340 

40.76 

5.1056 

1.0 

1.4938 

24.17 

4.497 

2.0 

2.7078 

14.34 

2.9075 

5.0 

4.2546 

7.7261 

2.4816 

10.0 

5.0487 

4.2665 

2.1617 

20.0 

5.3265 

2.1758 

1.9844 

100.0 

5.4211 

io.4373 

1.9146 
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TABLE  2 


Load,  Attitude  Angle  and  Threshold  Speed  for 

MC^ 


0.1 


0.2 


0.4 


0.6 


0.8 


0.9 


0.1 

0.00370 

88.29 

.2578 

.4976 

0.5 

0.01831 

81.52 

.5559 

.4979 

1.0 

0.03551 

73.43 

.7556 

.4985 

2.0 

0.06394 

59.38 

.9387 

.4996 

5.0 

0.10551 

34.32 

1.0340 

.4997 

10.0 

0.12180 

18.90 

.9774 

.4996 

20.0 

0.12757 

9.72 

.9240 

.4998 

50.0 

0.12899 

3.92 

.9023 

.5000 

100.0 

0.12914 

1.96 

.8953 

.5000 

0.1 

0.00758 

88.11 

.3671 

.4902 

0.5 

0.03747 

80.68 

.7974 

.4912 

1.0 

*0.07249 

71.97 

1.0778 

.4935 

2.0 

0.12998 

57.48 

1.3388 

.4979 

5.0 

0.21650 

33.08 

1.4487 

.4986 

10.0 

0.25291 

18.25 

1.3751 

.4982 

20.0 

0.26547 

9.39 

1.2945 

.4992 

50.0 

0.26935 

3.79 

1.2693 

.4999 

100.0 

0.26992 

1.89 

1.2593 

.5000 

0.1 

0.01660 

87.34 

.5703 

.4567 

0.5 

0.08201 

77.07 

1.2386 

.4595 

1.0 

0.15829 

65.86 

1.6354 

.4667 

2.0 

0.28302 

49.96 

1.9101 

.4840 

5.0 

0.48529 

28.38 

2.0071 

.4918 

10.0 

0.58882 

15.85 

1.9051 

.4917 

20.0 

0.62882 

8.19 

1.7935 

.4965 

50.0 

0.64164 

3.31 

1.7304 

.4993 

100.0 

0.64353 

1.65 

1.7237 

.4998 

0.1 

0.02861 

85.75 

.85156 

.3968 

0.5 

0.14364 

69.90 

1.8162 

.3978 

1.0 

0.28529 

54.90 

2.3057 

.4014 

2.0 

0.52462 

38.24 

2.5309 

.4180 

5.0 

0.92022 

21.43 

2.4221 

.4536 

10.0 

1.16019 

12.33 

2.2424 

.4713 

20.0 

1.26773 

6.46 

2.0759 

.4888 

50.0 

1.30417 

2.62 

1.9867 

.4980 

100.0 

1.30965 

1.31 

1.9790 

.4995 

0.1 

0.04635 

82.29 

1.3671 

.3230 

0.5 

0.25550 

56.43 

2.9426 

.2984 

1.0 

0.56879 

38.32 

3.7185 

.2715 

2.0 

1.12187 

23.92 

3.6096 

.2909 

5.0 

1.96448 

13.13 

3.0526 

.3498 

10.0 

2.52521 

7.93 

2.4862 

.4214 

20.0 

2.82464 

4.26 

2.1657 

.4712 

50.0 

2.93462 

1.74 

100.0 

2.95150 

0.87 

0.1 

0.060917 

78.12 

1.9140 

.2867 

0.5 

0.39428 

44.07 

5.2711 

.1846 

Table  2  (con'td. ) 


€ 

A 

P 

m 

a 

0 

% 

aM  S 

a.(f) 

0.9 

1.0 

0.98910 

26.83 

7.2323 

.1414 

.2.0 

2.04802 

15.71 

4.3436 

.2354 

5.0 

3.54985 

8.41 

3.3341 

.3091 

10.0 

4.55718 

5.21 

2.6243 

.3892 

20.0 

5.14859 

2.86 

— 

50.0 

5.37737 

1.18 

100.0 

5.41302 

0.59 

— 

[ 

• 

00 

• 

1 

TABLE  3 

I 

1  • 

Load, 

Attitude  Angle 

and  Threahold  Speed  for 
MC  ^ 

L 

D 

-%  '' 

1  e 

A 

P 

ct 

1 

m 

0 

F 

j  0.1 

0.1 

0.00119 

89.45 

0.1484 

Mote ;  There 

1 

0.5 

0.00592 

87.24 

0.3432 

la  no  A  ■  50.0 

1.0 

0.01180 

84.50 

0.4881 

for  L/D  ■  ^ 

2.0 

0.02329 

79.11 

0.6787 

! 

5.0 

0.05357 

64.51 

0.9916 

10.0 

0.08673 

46.67 

1.1233 

20.0 

0.11302 

28.08 

1.0794 

100.0 

0.12843 

6.11 

0.9698 

0.2 

0.1 

0.00246 

89.37 

0.2266 

0.5 

0.01227 

86.85 

0.5027 

1.0 

0.02444 

83.73 

0.7148 

2.0 

0.04809 

77.68 

0.9884 

5.0 

0.1095 

62.14 

1.4131 

10.0 

0.1770 

44.55 

1.5842 

20.0 

0.2333 

26.79 

1.5285 

! 

100.0 

0.2682 

5.83 

1.3553 

0.4 

0.1 

0. 00564 

89.00 

0.3672 

0.5 

0.02816 

85.03 

0.8090 

1.0 

0.05600 

80.18 

1.1440 

i 

2.0 

0.10960 

71.27 

1.5461 

5.0 

0.2443 

52.42 

2.0836 

10.0 

0.3930 

36.43 

2.2220 

- 

20.0 

0.5345 

22.02 

2.0918 

100.0 

0.6383 

.4.841 

1.8548 

0.6 

0.1 

0.01050 

88.12 

0.5859 

j 

0.5 

0.05259 

80.69 

1.2909 

1.0 
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1 
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AUis-Ctelaera  NaaufaeUalaR  CcaRiy 
Mllvaukse  1,  wlseoMla 

J.  V.  Pecker  Dlvtslce 
Aasricaa  Optical  Co^aap 
0109  Raws  Blvd. 

Plttsburkh  13,  PtaBsylvaalB 

Aasricaa  Society  of  Ltd>rloatloo 
tad  seers 

5  Bortb  Vsbasb  Avesue 
CblcaRO  2,  lUlaols 

Chalraaa 

Researcb  CoaaitUe  tm  LiArlcatloa 
Ybe  Aasrlcsa  Ooelcty  of  Mis^ai sal 
ttdseeri 

Ualtod  rnfliml  br  Coator 
505  Bast  kru  iU-eot 
Bev  York  17,  Bsv  York 

Nr.  Jaass  R.  Kerr,  PrvsUast 

lyeoalaR  OivtslM 

AYCO 

Btraiford,  Coaaoottr.ut 

keeaarek  Preelstca  Msnbanisas 

Oivlsloa  of  tko  Bardea  Corporatloi 

0  Old  Bevtowa  R«aR 

OaaRiaT,  CoaBMUeat 

Atlas  0.  L.  tOm,  Tlee  Prooldeat  • 


Jadvstrial  Park 

Port  WasklaRtoa,  PoaMsylvi 

m«a  Mvlelaa 
fhe  BMdla  Bsrparatli 
tii  Roowa  Avsaos 
Ptleo,  Rev  Yerk 
Altai  0*  Russell  T, 


Arasd  Rervlces  Tecbaical  lafceaatli 

ArliaRtca  Ball  ftatioa 
Arliactca  12,  Virciaia 

Office  of  Asslstaat  Director 
(Anv  Xsactorc) 

Oivtsisa  of  Reactor  Dsvslopasai 
U.  0.  AtesMs  Mora  CoHdssisa 
WaaRlaRtea  25,  0.  C. 

Atui  Mr.  Claraase  t.  MUor,  Jr. 

Mr.  f.  OrossaHk,  ChUf 
rnalaacriil  Dovolopant  Rraasb 
Isaetor  Dovelopaait  DtvlslM 

U.  B.  Atoads  MMrcr  Coaslsslca 
Vasblactoa  25,  0.  C. 

Readeuarters  Llbrsay 

V.  0*  Alcaic  Rasrey  Coadseibo 
WaahlBRtOB,  0.  C* 

Chief,  Oivlsloa  of  00lRd>fiR0 
Naritlw  AdalalstrailOT 
OAO  BuIJJIbr 
WashlaRtai  25,  B*  C, 

CrysRMls  liiln«r1  >M  laboratory 
■atloaal  Rwoau  of  Btoadards 

Bouldor,  Colorado 
Attat  Library 

Mr.  Im’oM  loooiOR 


MaiaiotroUcB 

1512  B,  Biseet,  i.  R. 
WashlMtoa  25,  0.  C. 

Nr.  Bdaiad  R.  RUsoo 

Chief,  tumriaotiom  A  Voar  1 
Lsvls  Rosaareh  Caaisr 


Adsdaisiratiia 
21000  Broohpark  Road 
OlovsUad,  Ohio 


1 


Nr,  NtNolpli  B.  Bspof 
Ooord  e#  Marshall  ipaso  PUdii  flsator 
2  RailssMU  Aoroaaoiles  aad  Rpaso 

1  Adalalstrailca 

1  Ouidaass  aad  Ccairol  Olrisica 

1  Oyrpo.RUbUiaer  Rraasb 

X  MtatevlUs.  Alobaai  If 


Bo0U  Avlailcs  CorporailBN 
Bssoareb  laheratoiisa  OtrUi« 
Bowurtsld,  MsMiaa 
Atui  Nr.  Ralph  e.  Ura<ai 

Mr.  C.  B.  AOoaa 
Fhyslcs  TacbaelatP  PoparUHt 
AoTB  B0fs  Divijua 
Iho  RoelJR  CtaBay 
Boattlo  »,  VashURtaa 

Iryaat  rhirklaR  Orlador  Caaatf 
U)  CUaisa  leiaii 
BprUdfloU,  Varaoat 
Atui  m.  RaaU  cua 

Cadi  lias  ftaos  CoMaay 

P.  0.  Baa  3006 

Detroit  5,  1 - Irr 

AtU:  0.  J.  Tkgrlor,  Projoct  Haa 


Dr.  Dovay  J.  BaaloU 
Mraotar  of  Peeiiopaait 
Carrlor  Raaaareh  a 
Carrlar  FarRsap 
OyrwiM,  Bav  York 


Chawa  Toa^  CaepcraAlaa 

P.  0.  Rea  >907 


AtUi  ».  B.  C.  BlvUck 

YUe  ProaUaai  (MulaeerlaR) 

CbiTNlor  Corporatloa 
Refeaaa  Oparattaao 
P.  0.  Raa  737 
Detroit  31,  Mleblaas 
Atui  0,  C.  V,  Raldor 

J*  0. 

PiMoaltUB  NuUsore 
3pO0  MaU  Rtroet 
Boosica  ff,  faaas 


RHiiaaJ  Rtroao  of  I 
RoalRer,  Coierode 


OvUae  Wrldit  CcrporaAloa 

iwnMaTlial  MvIoIcb 


rtm  < 


t>v»trc,«  Hein* 

L4W61A  ieul«V4fd 
AA«*:«t  k^,  C4lifor«i« 

AttAi  Ag^rt  ff.  Saltk 

Bp*elal  ProjMt  MM^lr—r 

ford  iMtnaMSt 
il-'.O  TboMca  At^mm 
Uttf  ItUad  Citjr  1,  Vtv  York 
Atto'  Jarrli 

Mr.  Mol/  Ecu 
ford  Motor  Cc^pcajr 
SBcioMnoc  Md  Burr 

f.  0.  Box  2053 
bvorbori),  MloMfu 

Vr.  John  K.  M«/or,  Jr. 
lk«4llUlUCB  iMUai 
AppUod  Sciooeo  B«p«rtMat 
Bcloatirio  lAborotnry 
ford  Motor  Co^ptay 
f.  0.  Box  8053 
Doorboro,  Klchlpa 

AlB*t»4reb  NMufMturUc  DlviclM 
Tb«  OwTott  Curporciloa 
'>^51  U.  Oopuivoda  Boulevard 
l/f  Ax«»X«f,  OaUfontl* 

Attii(  J*r*y  Oit^nt,  Bupervlfor 

N*«ami1oc1  Ub<  •Depn,  yj.if 

OvMraJ  Ainmlu  DlvUioa 
OoMral  Dy.^MlQi  Corporation 
f.  0.  Box  fjob 
8m  D1«^  12,  CaU/oraU 

AtUi  Mr.  f.  V.  Olx^aoB  1 

BMTlac  Md  UBtrloaat  Cantor 
OMsrat  ItaclMortac  Laboratory 
OMoral  lloatrla  Co^iatiy 
1  Klvar  Bond 
tkbaoaotady,  Bow  York 

Attoi  0.  H.  fox,  NiMcor  2 

Mr.  L.  V.  VlM 

Oaoaral  Blactrla  Coapaajr 

Aircraft  Aeeaaaory  furblM  DopartMit 

9Vi  WaaUm  Avoaua,  BLU.  3*7k 

Limn,  Maaaacbuaatta  1 


Bataarah  Uberaterlaa 
OoMral  Motera  CcrperatlM 
OoMral  Notara  faohalaal  c«t«r 
U  Mila  and  Mw^  Boada 
1  Varraa,  MteblfM 

Attat  B.  Belaad  >kkl,  , 

NaahMtaal  OavalaixMM  Bopt.  i 

A.  C.  Bpark  riuc  Mnalea 
1  OoMral  NBtera  Corporatloa 
Nllvaukaa  1,  viaaonaia 
Attni  AllM  ttnUaaa  1 

Valtar  Oarov 
Xaarfott  Uvlaloe 
1  Oaooral  frMlalon  laa. 

U50  MaBrUa  avobm 

Uttla  ralla,  law  Jaraay  1 

OruMM  AlMraft  tadMailac  Caap. 
Batbpaca,  Lane  lalaod,  law  York 
Attat  If.  David  V.  Oraic#  Jr. 

Wwhiiadfcl  Oaal»  SmUm 
1  tBjlMarlu  DopartMt  1 

■ydroaciftlai,  iMorparaiad 
BiadaU  IdBaal  Bead 
levard  Oeuitr 

Uw«l,  HarjWd  1 

ItU/BdtlaMa  luilBaiB  NMhlBM  Oerp. 
Baaaarah  Ubepatary 

Bab  Jaaa,  CaUfortla 
AtUi  Or.  V*  B.  XnacLotA  8 


if.  1.  B.  Barr,  V.  f. 

Laar  Banaa  Divlalaa 
Abba  Bo^ 

lljTla,  Ohio  1 

Or.  Calua  0.  Ooataal,  0/53*30 
Bide.  201,  flaat  2,  faU  Alu 
Lockbaad  Nlaallaa  A  Spaea  Co. 

P.  0.  Box  50t 

Siauiyviaa,  CAllfornla  1 


Or.  J.  0,  ,  _ 

UttAB  BfetAM,  Zia. 

5500  Cabobb  AvABua 
Weedlaad  tflUa,  Call/crmU 


Mr.  Dca  Maara 

Uitaa  liataaa 


Vaedlaad  UUa,  CaUfaraU 


Mr.  A. 

Aatra  Dtrlalaa 
Ite  MiriMwdt  earporBUM 
1«555  fatuar  nraat 
Tab  fift,  CdUfArBU 


If.  XMiaU  Pandaa 
Yiaa  PraAldABC  (k#') 
MAOaBMU  AlTAfaft  OArpATAUAI 
TAMart  Bt.  •  ft.  to<da 


Mi^al^  Airpark 
Baa  516 

ft.  laulA  3,  MiAotvi 


Or.  Bane  ftarmllabl 
WaAhiBiaal  TAabBelay  iBaarparAlad 
1  Barbart  Orlva 
Utbaa,  Vav  Tark 


2 


1 


If.  I.  L.  BImm 

Itrtraalaa 

A  Mvlalas  aP  lartbrap  Oarparatiaa 
500  Baat  Oraa^tbaepa  Avaaaa 
AMbaU,  CaliferalA  X 

Bartraaiaa 

A  OlvlaUi  ar  lartbrap  Oarparatiaa 

100  Maraa  Btraat 

Boiveed,  Muaaabuaatta 

Atui  Mr.  B.  l.  jaaiaana,  faa.  Aaat. 

Praaialaa  Praduata  riianawi  l 


Pratt  B  Wdtaar  Airaraft 
1  Divlalaa  of  (MC  •  CAHL 

P.  0.  Boa  dU 
Middlatiw,  pnaaiatlat 
Attai  Ubrarlaa  X 

Ubtary,  BlAf.  lO-B-5 
Badla  CarperaUaa  af  Aaaflaa 
^  Caadaa  8,  law  Jaraajr  X 

If.  Babart  B.  Blaflar 
Baatftdiaa 

iualaaalaf  fiMlvlaUB 

0633  fiMap  Aval 

3  Caaap  Park,  OalifarMU  B 


Mr.  J.  tf.  Larar 
nilaf,  MpMap^BartiAl  0< 
lABayfaU  Aar#  Mvlalaa 
8600  BUper  Ba*d 
WlaaaapoUa,  MLaaraota 


Nr.  Carl  P.  Oraaaaar,  Jr. 
Dlraater  af  Baaaarah 
Bav  iaapablra  Ball  iaarlBp,  laa. 
Patarberaup,  I  ~ 


Nra.  Aliaa  bard,  UbrartM 

krdaa  Divlalaa  af  Obltad  Alraraft  Oarp. 

lAlaa  Btraat 

Borvalk,  CcBBAcilAiit 

Bortbrep  Cerporatlea 

Borair  Dinaloa 

lOOl  Baat  Broadvay 

Bavtbona,  Calif enla 

Alta}  Taekalcal  lafomatloa,  3125 


Bpiw  Aaraaapiial  Oabnap 

AttBi  MipBffrlaj  Ubrvp 

Uadbarp  PtaU 

Baa  DUp  18,  CaUfaraU 

Mr.  M.  A.  TraalaPl 
fABiarapB  B  Pertar 
78  Mall  Btraat 
Bav  Tark  5,  Bav  York 

Jaak  B  Baiatf,  A  Otvialaa  of 
tba  Uaflar  OarpavaUaa 
1725  ^  Btraat,  BatU  505 
Haablapaa  6,  D.  C. 

Mr.  Paul  A.  PIU,  Tiaa  PraaldMt 
tbClaMi'lac  B  Baaaarah 
Belar  Alraraft  CapHMp 
2200  bMlfia  UpMV 
Baa  Maa»  12,  CaUfarala 

If.  V,  0.  Vine 


0*2 

Oraat  Mack,  Wmr  Ycrk 


Ifa.  floranca  TurabuU 
M«itMwrlBC  Ubrarlaa 
Oprry  Ofroaeopa  Ccapaar 
Oroat  Bach,  Bav  York 

OuPairaad  AvUtloa  Daavar 
2bao  Vaat  70tb  AvaaiM 
Daavar  81,  Colarada 

tf .  Jam  V.  lUlMil,  PraaldMt 
nrboeraft,  laa. 

1<>6  B.  b/^la  Avaaua 
Moaroala,  CaUfoavla 

Ubivaraal  Maiah  Carporatlaa 
Avlmlaa  Ovpartaaat  Taofaaal  Ubraip 
kkOI  Cook  Avaaua 
Ut.  Loula  13,  Nlaaotrl 

Vaukaaha  Baariaiia  Corporatloa 

P.  0.  Box  )b6 

Maukaaha,  Vlaaoaala 

AtUi  Mr.  /.  M,  (f<for,  Cb.  tap. 

Mr.  Joba  Boyd 

Moaiiapouaa  Blactrla  Carporatlaa 

Baaoareb  Laboratarloa 

iMl  nttaburg,  Paaaaylvaalo 

Mr.  k.  Valtar 
Dlraetar  of  Paaaarcb 
Vorthiactaa  Carporatlaa 
borrlaoa,  Bav  Aaraay 

Dr.  V.  A.  Oroaa 
Aaaia  Corporatiaa 
9Jb  Obartar  Ilvd. 

Bodvood  City,  OaUfarala 

Prafaaaar  J.  Modray 

DaparlMat  of  Maabaal  ral  "-gl — liaa 

Ualaa  Collap 

Hbaaaatgy  B,  Bav  York 

Btratoa  Divlalaa 
faireiaiA  Btratoa  Carporatlaa 
Baytaara,  L.  1.,  Bav  Yark 
AtUi  m,  .fpa  NMABar 

A.  0.  BpPt  PlMI  MHbIp 

OAMfbl  Hatara  OarpArallaa 

■avU  U$ 

VabafUU,  taaaaxbaaatU 
AiUt  Taabalaal  Ubrary 


Tba  Clavalaad  Orapblta  Broaia  Co^aisy 
ITOOO  Bt.  Clair  Avaaua 
Clavalaad  10.  Ohio 

1  AtUt  Mr.  B«  B.  Jeaapbaea  1 

Mr.  Biahard  J.  tali,  Nwapr 
iaarlBB  Davalopaaot  B  Caatraat 
1  Bav  Dvpartvra  Divlalaa 

Oaaaral  Natara  Cerporatlea 

•rlatol,  Caaaaatlout  1 

If.  Jaaaa  J.  Barkar 
1  10  ValAM  Avaaua 

Jarlaba,  Saw  York  1 

Tbeaaaa  Bane  Uoeldrldca 
TAPOO  Croup 

1  Bav  Davleaa  laboraterlae 
7209  Platt  AvMua 
Clavalaad  k,  Ohio 

Attai  Mr.  0.  Oackar  1 

1  Dr.  L.  Uabt 

lataraatleaal  Buaiaaaa  Maablava  Corp. 
Tbeaca  J.  Vataoa  naaaareh  Caatar 
P.  u.  Bea  218 

Yorktava  Balpta,  Saw  York  1 

Mr.  B.  A.  Boblay 
taloa  Caiblda  Sueloar  Co. 
feat  Offlea  Bea  P 

Oak  Bldp,  ToMaeaae  1 

1 

leeMiaad  Alraraft  Corporatiaa 

Nlaallaa  Md  Bpaaa  Mvlalaa 
Taabalaal  lafaraatlaa  Oaavar 
3851  laaavar  Btraat 

>  Pala  AlU,  CaUforaU  1 

Offlaa  af  foakalaal  Barvlaaa 

BaparUtat  af  CaaiiiPia 

Vaablaftaa  15,  D.  C.  1 


1 


